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ABSTRACT 


Growth rates and parameter dependences are calculated 
for the current-driven instabilities of electrostatic (with 
finite-/? corrections) and electromagnetic ion cyclotron waves. 
For 0.25 < T /T. < 2.5, ion cyclotron waves have large 

G X 

growth rates, while ion acoustic waves are still stable. In 
fusion devices, where electrostatic waves may be stable, 
electromagnetic ion cyclotron waves are unstable for 

_3 

/3^ > 10 . The suggestion is made that ion cyclotron waves 

may provide an anomalous resistance, with 


V.. ~ A 1//2 ft . (T ,,/T. .). 
ef f ci e , II i,l 
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INTRODUCTION 


Ion cyclotron instabilities, driven by magnetic-field-aligned currents 

and/or electron heat conduction, are of interest to the physics of the solar 

1 2 3-5 

wind, the ionosphere, Q-machine plasmas, and to certain controlled 

thermonuclear fusion devices such as tokamaks, which have confining electron 
currents. We report recent salient results concerning four basic extensions 
of the Drummond-Rosenbluth^ investigation of the electrostatic ion cyclotron 
instability: arbitrary T^/T. , where T a / T i is the electron -to-ion tempera- 
ture ratio; growth rates well above marginal stability; effects of contamina- 
tions by low charge -to -mass ratio ions; finite /3. corrections, where /3. is 
the ratio of ion particle to magnetic field pressure. Moreover, we report 

some preliminary results concerning the electromagnetic ion cyclotron cur- 
7 8 

rent instability, ’ and discuss a few implications of saturated ion cyclotron 

9 

turbulence. The details of this work will be published elsewhere. 

ELECTROSTATIC ION CYCLOTRON INSTABILITIES 


For purposes of analysis, we have chosen isotropic Maxwellian distribu- 
tions (albeit with different T^ and T. ) drifting along a uniform-background 

magnetic field B e . Thus, we have neglected the effect of the self-consis- 
o z 

tent magnetic field created by the current. To facilitate calculation in 

analytically difficult parameter regimes, a computer program has been 

developed which solves the dispersion relation for small-amplitude electro- 
G 7 

static waves. ’ The general electrostatic dispersion relation is cast into 
the Gordeyev integral form. ^ Complex roots of the dispersion relation are 
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11 


obtained numerically by a standard algorithm, such as Muller's method. 
A given wave mode is found either by making a suitably accurate guess or 
by following known wave modes into the parameter range of interest. 


SINGLE ION SPECIES 

2 £ g 

A marginal stability analysis ’ * indicates that when T^/T. = 1, the 

2 2 

wave that is marginally stable to the smallest current has K^R. ~ 1, 

OJ 1. 2 , and K . /K = 12, where R. is the thermal ion Larmor radius , 

C 1 II X « 

CO is the wave frequency, and K^/K^ is the ratio of wave numbers perpendic 

ular and parallel to B . Here, we focus on unstable ion cyclotron waves. 

6 7 

The general electrostatic dispersion relation * has dependences upon the 

2 2 

following parameters: K^R. , K^/K^, V^/a^ .fthe electron drift relative to 

ions in units (2T /m )^^], m/M, T /T., and the ratio of ion plasma to 


e’ e 


e i 


cyclotron frequency GO .• We can eliminate some deoendences: We 

pr ci 

2 2 2 2 ^ 

consider CO ./^ . > 10 . since this implies weak c 0 . /Q . dependence, 
pi ci- pi ci 

We eliminate the dependences upon K^/K^ and K^R^ by searching numeri- 
cally for local maxima in the growth rate. Choosing M/m = 183 6, corre- 
sponding to an H* plasma, we plot in Fig. 1 contours of maximum growth 
rate for the n = 1 wave as a function of the two remaining parameters, 
V^/a^ and T^/T.. For comparison, we also plot the marginal stability 

contour for the = 0 ion acoustic wave, corresponding to the magnetic 

12 13 

field-free calculations of Fried and Gould and Stringer, and the contour 

for y /w . = 1, which marks the transition to the nonresonant electron- 
max' pi 

ion beam instability originally studied by Buneman. The shaded region 
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denotes the portion of the T /T., VL/a parameter space for which the ion 

e 1 D e 

cyclotron wave is unstable but the ion acoustic wave is stable. The ion 

cyclotron wave is unstable to the smaller current over the range 

0.05 < T /T < 8.5. Furthermore, it has a significant growth rate 
e i 

y /f2. =0.1 (while the ion acoustic wave is still stable) over the range 

0.25 < T /T. < 2.5. Thus, the ion cyclotron instability dominates the ion 

acoustic instability over an interestingly large portion of parameter space. 

It can be shown that the region of ion cyclotron dominance increases as 
2 

M. increases. 
i 

CYCLOTRON INSTABILITIES IN MULTI-ION PLASMAS 

We now turn briefly to the operation of the electrostatic ion cyclotron 

instability in ion mixtures. Besides the above-listed parameters, there is 

now an additional dependence on N^/N^ the relative ion composition. Only 

the marginal stability (y /tt . = 0) contour is considered here because of 

rusix c i 

the added dimensionality in parameter space. We may ask the following 
questions: First, in a mixture, the ion harmonic of which ion is unstable to 
the smallest current? Second, how do critical electron drifts respond to 
changes in ion composition? Third, at what point does the plasma switch ion 
instabilities as the ion composition is varied? fourth, and last, how do these 

effects vary with T /T.? Figure 2 shows the results of an analytic calcula- 

2 + 
tion of the electron drift V for marginal stability in a mixture of H and 

c 

0 + (a case of ionospheric interest) as a function of T /T. and fractional O^" 

e i 

concentration, N /N . It is convenient to represent our results in N /N , 
o+ e o+ e 
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V^/a space with T /T. as a parameter. Solid and dashed lines denote 

the marginal stability curves for the instabilities near the first O and H 

cyclotron harmonics , respectively. The appropriate temperature ratio 

+ + 

T /T. is indicated at the intersections of the H and O curves. For 
e x 

simplicity, the H and O temperatures were taken equal. Figure 2 
indicates that the heavy-ion harmonic is unstable to the smaller current for 
surprisingly small heavy-ion contaminations, of the order of 10%. 

ELECTROMAGNETIC ION CYCLOTRON WAVES 
A computer program which solves the full electromagnetic wave disper- 
sion relation for spatially homogeneous, drifting electron and ion Maxwellian 
distributions in a uniform magnetic field has also been developed. In addi- 
tion to the single ion parameters listed above, the dispersion relation now 

7 

depends upon /3^. For large-enough /3^, it has been shown by Stix that elec- 
tromagnetic ion cyclotron waves (which have electrostatic corrections at 
oblique propagation) have a parallel phase velocity sufficiently slow that they 

can be destabilized by small electron drifts. Details of our investigation of 

9 

this instability will be considered elsewhere. 

Figure 3 shows some of our preliminary results. For M/m = 1836 and 

T /T. = 1. we plot contours of constant maximum growth rate in a /3. , 

6 1 ; 1 

V^/a^ parameter space. 

For purposes o.f comparison, we have also plotted the marginal stability 

contour, y /ft. = 0, for the "electrostatic" ion cyclotron wave. This 
'max' i 7 

mode is electrostatically polarized for j3. small, and develops small 
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electromagnetic components as /3. increases. These finite j3^ corrections 
are stabilizing. However, ion cyclotron waves of one sort or another con- 
tinue to be unstable as B. increases, because of the lower threshold of this 

1 

-3 

new electromagnetic instability for /3^ > 10 . It is probably this mode 

-2 

which is most important in the solar wind, since /3. > 10 , even near the 

sun. 

. _2 

For V max A\ = 10" t the real frequency is typically 0.75fL, and 

K.R. ~1.0 -1.5, K..R. < ,1 for maximum growth . As the drift velocity is 
1 i II l 

_ 3 

decreased to the point where y /Q,. < 10 , K..R. for maximum growth 

max' i I l 

drops rapidly below 1.0 and K^R. increases somewhat above 0.1. Drifts 
then tend toward Stix's estimate. 

Before closing the discussion of the electromagnetic current instability, 
let us note that a drifting electron Maxwellian certainly is not the self-con- 
sistent equilibrium distribution in a finite j3. current-carrying plasma. 
However, the Maxwellian should reasonably well describe the electron 
growth. Furthermore, any stabilizing effect of the inhomogeneous magnetic 
field, created by the direct current, might be unimportant because of the 
small wave group velocity. Kan has indicated to us that if the inhomogeneous 

magnetic field determines a pressure gradient, a drift mode might also be 

_2 

easily destabilized for B. > 10 

i 

NONLINEAR CONSIDERATIONS 

One -dimensional quasilinear electron plateau formation can saturate the 

_2 

electrostatic ion cyclotron instability at e^/T^ <10 , 0 being the wave 
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potential integrated over the unstable spectrum. Dum and Dupree, 


however, have shown that if electron plateau formation is inhibited, 

3-5,15 

e$/T - 1/2. This large saturation level agrees with Q machine results, 


where end effects may prevent plateau formation. We have observed 


16 


that largo e^»/T permits ion cyclotron waves to interact with a significant 

portion of the electron distribution, at an effective electron-ion frequency 

V .. — A£2 (T ,,/T. \ ), where A is the atomic number, and 1 and II 
eff ci e , II 1 , x 

refer to B . For large e<b / T , anomalous resistance due directly to wave- 
o : e 

electron scattering may occur, or turbulence may hold back the bulk of the 
electron distribution while high-energy electrons carry current demanded 
above cyclotron threshold. In the latter case, anomalous resistance ulti- 
mately may be derived from two-stream instabilities generated by the elec- 
tron tail. Significant spatial plasma diffusion and perpendicular ion heating 
may also occur. 

The conclusions described above could be important for tokamak con- 

trolled-fusion devices. Experiments in the Princeton Model ST may be 

operating in regimes where the electrostatic ion cyclotron wave is likely 

17 -19 

to be unstable, but where the acoustic wave may or may not be unstable. 

It has been suggested that larger /3. (higher density) may eliminate electro- 
static current instabilities by reducing the electron drift needed to produce 
the confining magnetic field. However, we must still contend with the 
electromagnetic ion cyclotron instability. Bers, Manheimer, and Coppi 

have already suggested that ion cyclotron waves may account for the anoraa- 

• 20 

lous perpendicular ion heating observed in the Soviet TM-3 Tokamak. 
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FIGURE CAPTIONS 


Fig. 1. Contours of maximum growth for the n = 1 H cyclotron wave 

(solid) and the H + acoustic wave (dashed). The shaded region denotes 

where the cyclotron wave is unstable but the acoustic wave stable. 

Fig. 2. Critical drift versus fractional O concentration for 

T /T. = 1,3,5,10. The solid (dashed) lines denotes the n = 1, 0 + (H + ) 

' 4 * 

critical drift. The transition from O to H waves shifts from 
N(0 + )/N = 0.15, at T /T. = 1, to N(0 + )/N = 0.06, at T /T. = 10. 

6 6 1 6 0 X 

Fig. 3. "Electrostatic" and electromagnetic ion cyclotron maximum 
growth rate contours for M/m = 1836 and T /T. = 1. The electrostatic mode 
has Re(o>/£2J> 1, whereas Re(co/f2J < 1 for the electromagnetic mode. The 
dashed extension occurs because of the difficulty of distinguishing the 

_3 

electrostatic from the electromagnetic modes. Above /3. - 10 , the elec- 

tromagnetic mode is unstable to the smaller current. 
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